DFT calculation of the intermolecular exchange interaction in the magnetic Mn 4 dimer 
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The dimeric form of the single-molecule magnet [Mn403CLi(02CEt)3(py)3]2 recently revealed 
interesting phenomena: no quantum tunneling at zero field and tunneling before magnetic field 
reversal. This is attributed to substantial antiferromagnetic exchange interaction between different 
monomers. The intermolecular exchange interaction, electronic structure and magnetic properties 
of this molecular magnet are calculated using density-functional theory within generalized-gradient 
approximation. Calculations are in good agreement with experiment. 

PACS numbers: 75.50.Xx, 75.45.+j, 75.30.Gw, 75.30.Et 
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Single-molecule magnets (SMMs), such as 
[Mn 12 Oi2(CH 3 COO)i 6 (H 2 0)4]-2(CH 3 COOH)-4(H 2 0) 
(hereafter Mm 2 )i and [Fe 8 2 (OH)i 2 (tacn) 6 ]Br 8 -9(H 2 0) 
(hereafter Fes)^ have received tremendous attention due 
to macroscopic quantum tunneling 3 and possible use as 
nanomagnetic storage devices. Hysteresis loop measure- 
ments on the SMMs Mn 12 and Fe§ showed magnetization 
steps at low temperatures upon magnetic field reversal^ 
This is due to quantum tunneling between spin-up states 
and spin-down states despite a large effective spin 5=10 
for each molecule. The resonant tunneling fields in these 
systems are primarily determined by the magnetomolec- 
ular anisotropy. Recently a dimerized single-molecule 
magnet [Mn 4 3 Cl4(0 2 CEt)3(py)3] 2 (hereafter Mn 4 
dimer) where Et=CH 2 CH3 and py=NCsH5, has been 
formed 5,6 which exhibited qualitatively different tunnel- 
ing behavior: quantum tunneling prior to magnetic field 
reversal and an absence of quantum tunneling at zero 
field in contrast to other SMMs such as Mni 2 and Fe 8 pi 
To understand the basis for the qualitative deviation we 
have calculated both the magnetomolecular anisotropy 
and the intermolecular exchange interaction in the Mn 4 
dimer using density-functional theory. Our results con- 
firm that there exists an appreciable antiferromagnetic 
exchange interaction between monomers and that tun- 
neling fields in this dimer are strongly influenced by the 
presence of the monomer-monomer exchange interaction. 
This interaction produces a bias field that encourages 
monomeric magnetic-moment reversal below zero field 
and prevents two monomers from simultaneously flipping 
their magnetic moments at zero field. We determine that 
the origin of the exchange interaction is not dominated 
by either kinetic or exchange-correlation terms and that 
the total "exchange" interaction is in fact an order of 
magnitude smaller than the kinetic contribution. For 
Mni 2 and Fes, the intermolecular exchange interaction 
has not been observed experimentally and it is generally 
accepted that the overlap between neighboring molecules 
is negligible. 

In this work, we discuss calculations on the Mn 4 dimer 
which is formed by inversion of the three-fold symmet- 




FIG. 1: Mn4 dimer geometry. The dimer is formed by in- 
version of the threefold symmetric monomer. Each monomer 
has a magnetic core consisting of three ferromagnetically cou- 
pled Mn 3+ spins (S=2) coupled antiferromagnetically to one 
Mn 4+ spin (5=3/2) ion leading to a total spin of S = 9/2. 
The distance between the two central CI atoms marked as the 
dotted line was measured to be 3.86 A. 



ric monomer shown in FigQ The magnetic core of 
the Mn 4 monomer consists of three ferromagnetically 
coupled Mn 3+ (5=2) ions coupled antiferromagnetically 
to the remaining Mn 4+ (5=3/2) ion leading to a to- 
tal ground-state spin of 5=2x3—3/2x1=9/2 (refer to 
FigHJ. The core has a similar cubane structure as the 
inner core of the SMM Mni 2 , although there are four 
Mn 4+ for Mni 2 . We investigate the electronic structure 
and magnetic properties of this SMM Mn 4 using density- 
functional theory (DFT). We calculate optimized geome- 
tries for the Mn 4 monomer and dimer, their binding en- 
ergy, the monomeric magnetic anisotropy barrier (MAE), 
and the exchange coupling constant between monomers. 
Results are compared with experiment. 
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Our DFT calculations 7 are performed with the all- 
electron Gaussian-orbital-based Naval Research Labo- 
ratory Molecular Orbital Library (NRLMOL)£ Here 
we use the Perdew-Burke-Ernzerhof (PBE) generalized- 
gradient approximation (GGA). 9 Before discussing ener- 
getics and magnetic phenomena we discuss two structural 
issues. First, we have considered monomers and dimcrs 
that are terminated by both H and by the CH2CH3 
radicals found in the experimental structure. Second 
we have considered structures based on two conform- 
ers of the monomeric unit. While a complete vibra- 
tional analysis will be discussed in a later publication, 
all indications are that both conformers are stable. The 
conformers have slightly different arrangements of the 
pyridine ligands. The first conformer was identified by 
a density-functional-based geometry optimization of the 
hydrogenated monomer. The second conformer was iden- 
tified by improvements on the monomer deduced from 
the experimental x-ray data. In the remainder of the pa- 
per we refer to these monomers as the computationally 
determined conformer (CDC) and the experimentally de- 
termined conformer (EDC). 

Each Mn4 monomer has threefold symmetry so there 
are 26 inequivalent atoms to consider. The number of 
inequivalcnt atoms is reduced to 20 when the CH2CH3 
radical is replaced by H. A pyridine ring is initially con- 
structed to lie in the plane defined by the vector con- 
necting Mn 3+ and neighboring N and the sum of the 
two vectors connecting Mn 3+ with two closest CPs (re- 
fer to Fig. ^) . The geometries for the pyridine ring and 
the cubane were first optimized separately to generate 
an initial geometry for the DFT calculations on the full 
monomer. The initial geometry for the monomer was 
relaxed using NRLMOL with the CI atom fixed to re- 
produce the experimental Cl-Cl distance (3.86 A) upon 
dimerization (i.e. adding inversion symmetry). Relax- 
ation continues until forces exerted on all atoms become 
^0.001 hartree/bohr. The CDC dimer is then obtained 
by inversion of the CDC monomer with the fixed value 
of d = 3.86 A (marked as dotted in Fig. Q). For the case 
of the x-ray deduced experimental geometry, the C-H 
bond lengths are underestimated (0.71 A to 0.96 A) in 
comparison to standard hydrogen bond lengths, which 
yields self-consistent forces on hydrogen atoms as large 
as 0.8 hartree/bohr. To improve the experimental geome- 
try, all hydrogen positions were first moved to create C-H 
bond lengths as 1.1 A, and then additional optimization 
of the experimental geometry was performed with the 
fixed Cl-Cl distance. The experimental geometry without 
corrected hydrogen positions was 53 eV higher in energy 
than that of the structure with corrected hydrogen po- 
sitions. Hereafter unless we specify, the EDC monomer 
refers to the optimized experimental geometry with cor- 
rected hydrogen positions. 

We have used full basis sets for all six different atoms 
and fine mesh. 10 Charges and magnetic moments for 
Mn's from the CDC monomer agree well with those from 
the EDC monomer. For example, a sphere with a ra- 
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FIG. 2: Electronic density of states (DOS) for majority and 
minority spins for the Mn4 monomer. Shown are projected 
Mn(3d) DOS of the two types of Mn ions, projected p DOS 
of the three N atoms and the four CI atoms, projected p 
DOS of the nine O atoms, and the total DOS defined by 
the sum of projected DOS of all atoms in the monomer. All 
projected DOS have the same scale which is different from 
that for the total DOS. The vertical line denotes the Fermi 
level. Directly below the Fermi level, for majority spins the 
projected Mn 3+ (3d) DOS has more weis ht than the Mn 4+ (3d) 
DOS. For minority spins, the tendency is the opposite. 

dius of 2.23 Bohr captures charges of 23.4 and 23.7, and 
magnetic moments of 3.6/is and — 2.5/is for Mn 3+ and 
Mn 4+ respectively. The total magnetic moment for the 
monomer is 9/is in good agreement with experiment. 
The HOMO-LUMO gap for majority (minority) spin is 
1.02 eV (2.42 eV). The energy difference between the 
minority (majority) LUMO and the majority (minority) 
HOMO is 1.17 eV (2.28 eV), which ensures that the sys- 
tem is stable with respect to the total magnetic moment. 
As clearly seen in Fig. [21 right below the Fermi level for 
majority spins the projected Mn 3+ (3cQ DOS is dominant 
over the projected Mn 4+ (3d) DOS, while for minority 
spin the opposite trend is observed. This confirms the 
experimental picture of three Mn 3+ spins antiferromag- 
netically coupled to a Mn 4+ spin. 

We calculate the binding energy by subtracting the 
dimer energy from twice the monomer energy. We find 
that the dimer is stable for both the CDC and EDC. For 
the CDC (EDC), the binding energy is about 0.16 eV 
(0.78 eV). The magnitude of the binding energy sug- 
gests attractive electrostatic interactions between differ- 
ent monomers. The discrepancy between the binding 
energy for the CDC and that for the EDC may be at- 
tributed to our substitution of ethyl for hydrogen in the 
CDC and/or the fact that the plane where a pyridine 
ring sits is different for both geometries. To check the 
former possibility, we calculate the binding energy of the 
EDC terminated by hydrogen, and obtain 0.45 eV. We 
have also verified that the conformation of a pyridine ring 
for the EDC is slightly different from that for the CDC. 
Thus, the discrepancy arises from both reasons. 

We have calculated the monomeric MAE in zero mag- 
netic field for both the CDC and EDC with the assump- 
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TABLE I: Binding energy, monomeric magnetic anisotropy 
barrier (MAE), and antiferromagnetic exchange constant J 
for the CDC with the distance between the two central 
Cl's held as the experimental value, d = 3.86 A [DFT(l)], 
the EDC with d = 3.86 A [DFT(2)], and the same as 
DFT(2) except that ethyl is replaced by hydrogen [DFT(3)]. 
DFT(4), DFT(5), and DFT(6) denote the same as DFT(2) 
except that d=3.86 A+l Bohr, d=3.86 A-0.5 Bohr, and 
d=3.86 A— 1 Bohr, respectively. The experimental values are 
from Ref@. The numerical uncertainty in the estimated val- 
ues of J is ~ 0.04 K. 
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tion that spin-orbit coupling is a major contribution to 
the MAE. For this calculation, we follow the procedure 
developed in RefUll Our calculations show that the Mn4 
monomer has uniaxial anisotropy along the threefold axis 
(the bond between Mn 4+ and CI in the cubane) , in agree- 
ment with experiment For uniaxial systems, the en- 
ergy shift A due to the spin-orbit interaction can be 
simplified to —^ ZZ {S Z ) 2 up to constant terms indepen- 
dent of (S z ) if the z axis is assigned as the easy axis»ii 
Then the classical barrier (MAE) to be overcome to 
monomer magnetization reversal M z —+9/2 to M z =— 9/2 
is 7 Z2 ((9/2) 2 - (1/2) 2 ). For the CDC (EDC) monomer, 
the MAE is 11.3 K (11.6 K), which is close to that for 
the hydrogenated EDC monomer. As shown in Table [I] 
all these numbers are close to the experimental value of 
14.4 K. The difference between our estimated MAE and 
the experimental value might be ascribed to other effects 
on the barrier such as spin-vibron coupling^ 

To calculate the exchange coupling constant J be- 
tween monomers, we assume that a monomer is an ideal 
5 = 9/2 object and that its effective spin is aligned 
along the easy axis and of Ising type (either M 2 =+9/2 
or —9/2). Then we calculate self-consistently energies 
of ferromagnetic (parallel monomeric spins) and antifer- 
romagnetic configuration (antiparallel monomeric spins) 
of the dimer, and take a difference 8 between the two 
energies. We find that the antiferromagnetic configura- 
tion is favored. The antiferromagnetic exchange constant 
J is determined from S = 2J(9/2) 2 . For the CDC, the 
energy difference is 31 microhartree so that J=0.24 K, 
while for the EDC, J=0.27 K. These can be compared to 
the experimentally measured value of J— 0.1 K£ The nu- 
merical uncertainty in the total-energy difference for our 
DFT calculations is at most 5 microhartree, which can be 
translated to the uncertainty in the exchange J as 0.04 K. 




[CI-CI distance] - exp value (3.86 A) 



FIG. 3: Logarithm of exchange constant J as a function of 
the monomer-monomer distance relative to the experimental 
value. The numerical uncertainty in J is ~ 0.04 K. The slope 
of the curve is about -2. 



We achieve high-accuracy in the total-energy difference, 
because we use exactly the same optimized dimer geom- 
etry with the same parameter values for a self-consistent 
approximation except for the effective spin configura- 
tions of monomers. Although our DFT estimated value 
of J is somewhat higher than the experimental value, 
this may be acceptable considering the assumptions we 
made and the fact that DFT calculations often overes- 
timate exchange interactions. In some cases, the PBE 
generalized-gradient approximation may not fully cancel 
the self-interaction in the Coulomb potential. Therefore, 
the electrons in our calculations are slightly more diffuse, 
which should lead to overestimated exchange interaction. 

It is interesting to examine whether the exchange in- 
teraction varies significantly with the monomer- monomer 
separation. We consider the case that each monomer 
is displaced toward or away from the center of mass of 
the dimer along the easy axis. Then we calculate the 
exchange constant J for the EDC dimer with three dif- 
ferent monomer-monomer distances from the experimen- 
tally measured value. The monomer-monomer distance 
is varied by changing the two central Cl-Cl distance with 
a monomer geometry fixed. If the central Cl-Cl bond 
length increases by 1 Bohr, then J decreases down to 
0.10 K. If the bond length decreases by 0.5 Bohr (1 Bohr), 
J increases to 0.47 K (0.81 K). Table [I] summarizes the 
separation dependence of J and of the monomeric MAE. 
As shown in Table [I] the monomeric MAE does not de- 
pend on the exchange interaction between monomers, be- 
cause the monomer geometry has not changed during this 
process. Figure |21 shows that J increases exponentially 
with decreasing the separation distance. This tells us 
how quickly the overlaps of neighboring wavefunctions 
decrease with increasing the distance. We have decom- 
posed the J values into kinetic, coulombic and exchange- 
correlation contributions. The kinetic contribution is an 
order of magnitude larger than the total value of J and 
it is significantly cancelled by the exchange-correlation 
contributions to the J value. 

Since we estimated the anisotropy barrier and ex- 
change constant, we can construct a model Hamiltonian 
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TABLE II: Initial \M 1 ,M 2 ) and final states \M[,M 2 ) partic- 
ipating in quantum tunneling at resonant fields, B vca . This 
was calculated by exact diagonalization of Hamiltonian Q. 
Mi and M 2 are the projected magnetic moment along the 
easy axis for each monomer. For clarity, degeneracy in ini- 
tial/final states is not listed. Only for the case with 7 ZZ 2> J 
the initial/final states are eigenstates of Hamiltonian Q. The 
third and fourth resonances are from one degenerate state to 
another degenerate state, and they are split due to transverse 
terms in the exchange interaction. The same logic is applied 
to the last two resonances. B^,* p is the resonant field for 
7 ZZ = 0.72 K and J = 0.1 K. B?* T is for 7 ZZ = 0.58 K and 
J = 0.27 K B^° d is for 7zz = 0.58 K and J = 0.1 K. 
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for the dimer according to 

n = -j zz (s 2 lz + si) + j ^ ■ s 2 (i) 

where the uniaxial anisotropy parameter j zz = 0.58 K 
and J — 0.27 K. To determine whether our values of 7 ZZ 
and J can reproduce the experimental values of the reso- 
nant tunneling fields (Fig. 4 in Ref@), we calculate these 
fields using exact diagonalization of Hamiltonian . Al- 
though it is crucial to include some small transverse 
terms in the Hamiltonian (such as transverse anisotropy 
and transverse fields) for calculations of tunnel splittings, 



the transverse terms do not affect the resonant fields 
much. Table ITT1 summarizes the resonant fields for some 
low-energy states for three different values of ^ zz and J: 

(1) the experimental values ^ zz — 0.72 K, J — 0.1 K; 

(2) 7 Z2 = 0.58 K, J = 0.27 K; (3) slightly modified 
version of our DFT results 7 ZZ = 0.58 K, J = 0.1 K. 
Let us focus on two tunnelings which were prominent 
in the experimental measurements: \M\ = 9/2, Mi = 
9/2) -> \Mi = 9/2, M 2 = -9/2), and 9/2, —9/2) -> 

— 9/2,-9/2), where Mi and M2 are the eigenvalues of 
the spin operator projected along the easy axis for each 
monomer. For these two tunnelings, the resonant fields 
are solely determined by J and are independent of ~f zz : 
B ICS ~ T9J/(2(7/ib). Therefore, model Hamiltonian 
with our estimated values will not quantitatively repro- 
duce the experimental resonant fields. However, in this 
case (when ^ zz becomes comparable to J), we notice that 
the hysteresis loop exhibits richer features such as more 
magnetization steps before magnetic field reversal. Since 
DFT often overestimates exchange interactions, we also 
calculate the resonant fields with J decreased to 0.1 K 
and 7 22 fixed to examine if agreement with experiment 
improves. We find that some resonances agree with ex- 
periment and some do not agree. 

In summary, we have calculated optimized geometries 
for a monomer and dimer of the SMM Mn4 using DFT. 
For both the CDC and EDC, we calculated binding en- 
ergy, monomcric MAE, and the exchange interaction 
between monomers. The binding interaction between 
monomers is electrostatic. Our calculated anisotropy 
barrier is close to the experimental value. The exchange 
interaction between monomers is twice or three times 
larger than the experimental value. Overall, our DFT 
calculations are in qualitative accord with experiment. 
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